Higgs-lepton inflation in the supersymmetric minimal seesaw model 
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We investigate a scenario of cosmological inflation realised along a flat direction of the minimal 
seesaw model embedded in supergravity with a noncanonical R-parity violating Kahler potential. 
It is shown that with appropriate seesaw parameters the model is consistent with the present ob- 
servation of the cosmological microwave background (CMB) as well as with the neutrino oscillation 
data. It is also shown that the baryon asymmetry of the Universe can be generated through leptoge- 
nesis. The model favours supersymmetry breaking with the gravitino as the lightest superparticle, 
and thus indicates the gravitino dark matter scenario. An interesting feature of this model is that 
the seesaw parameters are constrained by the CMB spectra. The 2-ct constraints from the 9-year 
WMAP data yield a mild lower bound on the seesaw mass scale > TeV. We expect that the obser- 
vation by the Planck satellite will soon provide more stringent constraints. The phenomenological 
and cosmological implications of the R-parity violation are also discussed. 
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I. INTRODUCTION 

In cosmology, the particle physics origin of inflation re- 
mains as an unsolved problem. In search of a realistic in- 
flationary model one may take two possible approaches. 
One is the top-down approach, starting with a unify- 
ing theory including gravity. Sring cosmology loop 
quantum cosmology , and cosmological model building 
based on F-theory grand unification [3] fall into this cate- 
gory. The top-down approach is an ambitious program to 
build a consistent scenario from the first principles. The 
other is the bottom-up approach, that is, to investigate 
embedding of inflation into particle physics models that 
are confirmed at low energies, starting from the Standard 
Model (SM) and its various extensions. The bottom- 
up approach is advantageous in predictability and falsi- 
fiability. In view of the remarkable progress and ever- 
increasing precision in observational cosmology 0, Q one 
could hope that in the near future observation will be 
able to narrow SM-based inflationary models down to a 
handful of candidates. In this paper we take the bottom- 
up approach and propose a model of inflation, which we 
believe to be a promising and testable candidate. 

Recently, there has been a revival of interest in cos- 
mological inflation within the SM, where the Higgs field, 
nonminimally coupled to gravity, is identified as an infia- 
ton LGrilSj (see also j20|J24'| for older proposals of nonmim- 
ilally coupled inflation models). Interestingly, it requires 
the Higgs mass to be 126 — 194 GeV, consistent with the 
mass scale suggested by the Large Hadron Collider exper- 
iments ^5j, i2(|| . The Higgs inflation model also predicts 
small tensor-mode perturbation that flts remarkably well 
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with the present-day cosmological microwave background 
(CMB) observation. The success of the model comes at 
the cost of an extremely large coupling ^ ~ 10'* between 
the Higgs fleld and the background curvature, which 
could lead to the violation of the unitarity bound [l^ - 
HH (for the controversy see also [iMl Hi HI]). Also, 
the model suffers from the usual hierarchy problem of the 
SM, namely the Higgs potential is destabilised due to ra- 
diative corrections^, |^, [H, [l^. As is well known, the 
hierarchy problem is solved, or at least tamed, by intro- 
ducing supersymmetry. Supersymmetric versions of the 
Higgs inflation model were considered in [27l - [29j . and it 
was shown that while embedding into the minimal super- 
symmetric Standard Model (MSSM) is not possible (as 
the fleld content is too restrictive), the next-to- minimal 
supersymmetric SM successfully accommodates Higgs in- 
flation. Embedding into supersymmetric grand unified 
theory [s^l (see also [3l| ) and other supersymmetric SMs 
[3^ have also been shown to be possible. 

Strictly speaking, the SM is not entirely a satisfactory 
model of particle phenomenology any longer, even apart 
from the hierarchy problem: it fails to explain neutrino 
oscillations; it does not contain good candidates for the 
dark matter; the generation of the baryon asymmetry in 
the electroweak phase transition is known to be prob- 
lematic. Among the simplest extensions of the SM that 
can solve all these problems is the supersymmetric exten- 
sion of the seesaw model (S^l • In [s^ we proposed a new 
scenario of inflation within the supersymmetric seesaw 
model, inspired by the developments of the nonminimally 
coupled Higgs inflation models. Notably, in that scenario 
the problem associated with the large nonminimal cou- 
pling is alleviated. The purpose of the present paper is 
to discuss the viability of this inflationary scenario using 
realistic neutrino mass parameters. We follow the mini- 
malistic guiding principle of Hig gs inflation and consider 
the minimal seesaw model [35|, i.e. with two families 
of the right-handed neutrinos. There are several other 
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inflationary scenarios based on the supersymmetric see- 
saw model [36r.42.] . We shall compare the cosmological 
predictions of these models with ours, and argue that at 
least some of these models are soon to be excluded by 
CMB observations. 

The main question we shall address below is whether or 
not the model allows parameter regions that are consis- 
tent with the observed baryon asymmetry, the spectrum 
of the CMB and the neutrino oscillation data. The ra- 
tio of the baryon density ub to the entropy density s is 
measured to be d, @ 

= — = 8.7 X 10-". (1) 

s 

As mentioned above, the electroweak baryogenesis within 
the SM is problematic as it requires rather implau- 
sible strongly first order electroweak phase transition. 
We show that in our scenario the thermal leptogenesis 
[43I with the reheating temperature Trh ~ 10^ GeV, 
combined with the resonant enhancement effects (43 - l46| 
yields the observed amount of the baryon asymmetry llj . 
An interesting feature of our scenario is that the CMB 
spectrum is related to the seesaw mass scale [sj]. We 
investigate this relation in the supersymmetric minimal 
seesaw model, and show that for a wide range of param- 
eters our scenario is consistent with the present CMB 
observation. In the near future the CMB data will fur- 
ther constrain the seesaw mass scale. In our scenario 
the R-parity needs to be broken. We discuss that this 
requirement leads to a scenario of gravitino dark matter. 

The paper is organised as follows. In the next section 
we describe the supersymmetric minimal seesaw model 
on which our scenario is based. In Section III the dy- 
namics of inflation is discussed, and in Section IV the 
reheating process and baryogenesis are studied. The re- 
lation between the cosmological parameters and the neu- 
trino mass parameters is discussed in Section V. Section 
VI deals with the R-parity violation. We conclude in Sec- 
tion VII with comments. Technicalities in computing the 
baryon asymmetry are relegated to the Appendix. 

II. THE SUPERSYMMETRIC MINIMAL 
SEESAW MODEL 

Our model is based on the supersymmetric seesaw 
model, with the superpotential, 

W = tiHuHa + yl^u'rQ.jH^+y'jd'iQjHd 



Here, Q, u'^, (P, L, e'^, Hd are the MSSM superfields, 
the right-handed neutrino superfields (having odd 
_R-parity), Mm the corresponding right-handed neutrino 
mass parameters, the MSSM /z-parameter and j/^-* , y'^ , 
2/^-', yl^ are the Yukawa couphngs. The superfields Q, L, 
Hu, Hd are SU{2) doublets and the contraction using 
the SU{2) invariant i(T2 = { °i i ) is implicit, whereas 
M^, d'^, e'^ are SU(2) singlets. In this paper we shall 
focus on the simplest realistic seesaw model with two 
right-handed neutrinos (the minimal seesaw model 35] ) . 
Thus the family indices run m, n, • • • = 1, 2 for the right- 
handed netrinos and i, j, • • • = 1, 2, 3 for the other lepton 
and the quark superfields. 

Using the stationarity condition SW/SN^ — 
y^'^LiHu + MjnN^ = 0, the superpotential along 
the flat direction reads 

= -lM,-'{yWL..H^){y7L,H^). (3) 

The Higgs field develops the vacuum expectation value 
at low energies. 



where (7J°) = ^sin^ with v = 246 GeV. We use 
tan /? = 10 throughout this paper. The mass matrix of 
the left-handed neutrinos obtained from ^ is 



'm-u = m^M ^rriD, (5) 



where m/) = yniH^^) and M = diag(Mi, M2). This is the 
celebrated seesaw relation. The Maki-Nakagawa-Sakata 
(MNS) lepton flavour mixing matrix is parametrized as 



C12C13 S12C13 sue''^ \ / 1 

t^MNS = I -S12C23 - Ci2S23Si3e"' C12C23 - Si2S23Sl3e"^ S23C13 6*'^ | , (6) 

S12S23 - Ci2C23Sl3e'^ -C12S23 - Sl2C23Sl3e"' C23C13 / \ 1 
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where Sij = sinOij, Cij = cos 9ij, and S, a are the CP- 
violating Dirac and Majorana phases. The neutrino mass 
matrix is diagonahsed by the MNS matrix as 



(7) 



where = diag(TOi, m2, 7713). We use the neutrino mass 
parameters from the oscillation data [13, 13 : 

sin^ 2012 = 0.87, sin^ 2023 = 1-0, sin^ 26*13 = 0.092, 
Am?2 = ri4 -m\^ 7.59 x 10"^eV^ 
Am^a = \ml -m\\^ 2.43 x lO^^eV^ 

In terms of these the neutrino masses are 



(8) 



mi = 0, m2 = Y Am^2, TO3 = y ^"^12 + ^"^23' (9) 
for the normal mass hierarchy (NH) and 



mi 



Am23 — Amfji ™2 = y/ Amjg, m3 — 0, (10) 



for the inverted mass hierarchy (IH). The neutrino mass 
matrix can be conveniently parametrized as p9l [50l | 



MNS' 



(11) 



where \J M = diag(vMi, \fM-i) and 








J1^2 



In ([TT|) . is a 2 X 2 orthogonal matrix 



(12) 



cosw smw 
— sinw cosw 



w G C. 



(13) 



We shall write the real and imaginary parts of w as 
= a + z6, a, 6 G M. A merit of the minimal seesaw 
model is its strong predictive power. Note that the mass 
matrix m/j (or equivalently the Dirac Yukawa coupling 
?/£)) contains 9 real degrees of freedom (6 complex minus 
3 phases), 5 of which are constrained by the oscillation 
data, namely the two neutrino masses and the three an- 
gles in ([5]). In addition, there are one Dirac and one Ma- 
jorana phases, and the remaining two degrees of freedom 
correspond to the choice of a and h. The right-handed 
neutrino masses Mi, M2 and the Dirac and Majorana 
phases are not fixed by the present experiments. In our 
scenario these parameters are subject to the constraints 
from the CMB, as discussed below. 



III. INFLATIONARY DYNAMICS AND 
COSMOLOGICAL PARAMETERS 

In this section we describe the construction of the in- 
flation model and discuss its prediction on cosmological 



parameters. The model, which is a multi-family general- 
isation of the one introduced in [sj] , has some similarity 
to the supersymmetric Higgs inflation models p7l - [30t . A 
notable difference from these models is that the nonmin- 
imal coupling ^ is allowed to take small values and the 
unitarity violation problem is thus alleviated. 



A. The model of inflation 

We assume that inflation takes place along one of the 
L-Hu D-flat directions and as an initial condition one of 
the three (scalar components of the) Li flelds, call it 
with k fixed, has a large expectation value and dominates 
the inflationary dynamics over the other L^'s. The D-flat 
direction along L^-Hu can be parametrized as 



V2\0 



""""Til" 



(14) 



with k fixed. Disregarding Q, u'^, d'^, e'^, that play no 
role during inflation, the superpotential ([2]) reads 



For the Kahler potential K 
noncanonical form, 



(15) 

3$ we choose a slightly 



\H^\' + m\^ + \N^\') 
c 



* = l-l{\Lk\ 

= 1- J(i^i' + i^ri' + i^2i') 

+ l{^^+^^) + ^-{\N^\* + \N.I\^). (16) 

The terms proportional to 7 yield nonminimal coupling 
between the inflaton and the background curvature, and 
the other noncanonical terms proportional to C have been 
introduced for controlling the inflaton trajectory. The su- 
pergravity scalar potential (in the Jordan frame) is com- 
puted from 113 and (HH) as (see (51|]). 



3Wi 



inf 



(17) 



where the repeated rt's are summed over, k is flxcd (no 
sum), and 



A. 



A, = 



E ■ 

m=l,2 



m=l,2 

E ■ 

m=l,2 



i-4C|iv;^.l' ' 

iy^^^2 _^ M^^c^) N^^ (1 - 2C m 



l-4C|iV;^ 



cm 



l-AC\N- 



(18) 
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The Dirac Yukawa couphng is 
V2 



V sin/3 



(19) 



and in p7)) the fields cp and iV^ are understood to be 
the scalar components. In the Einstein frame the scalar 
potential is 



V -^^ 



B. The inflaton trajectory 



(20) 



the slow roll parameters becomes 0(1)) before the com- 
plication starts. In the example considered here the slow 
roll terminates a.t tp = 0.0371, where and are still 
stabilised close to — 0. Fig. O shows the shape of 
the potential as a function of (p and ReiVf , where ImTVf , 
ReA^I, ImA^I are taken to be along the trajectory of Fig. 
[T] In this case^ the primordial tilt and the tensor-to- 
scalar ratio of the CMB are computed to be Ug = 0.968, 
r — 0.00296. The red spectrum with very small r is typ- 
ical of the nonminimally coupled Higgs inflation. The 
horizon exit of the comoving CMB scale takes place at 
(fi = 0.318. In computing the CMB spectrum the role 
played by C is minor, as ( has no effects once the inflaton 
trajectory is stabihsed along « 0. 



The model of inflation we consider is a system of three 
complex scalar fields cp, and iV^f . While the dynamics 
could consequently be quite involved in general, it turns 
out that under mild assumptions the model reduces to 
that of single-field slow roll infiation. This is by virtue of 
the noncanonical terms in the Kahler potential propor- 
tional to 

When Mm are large our model is similar to the su- 
persymmetric versions of SM Higgs inflation discussed 
in |27H3G |. Let us consider, as an example, the see- 
saw masses Mi = M2 = 10^^ GeV and choose for con- 
creteness the normal mass hierarchy of the neutrinos and 
a = 0, 6 = 1, (5 = 0, (T = 0. Then the Yukawa coupling 
are found to be 

y}j = 0.0465 + 0.0233i, yj^ = 0.0306 - 0.0354i, 

y]^ = 0.0434 + 0.106i, y]^ = 0.139 - 0.0331i, 
yi3 = -0.0536 + 0.106i, yj^ = 0.139 + 0.0408i. (21) 

The dynamics of inflation can be studied by examining 
the steepest descent trajectory of the scalar potential 
Assuming that inflation takes place in the Lk-Hu 
D-flat direction along the first generation lepton super- 
multiplet {k = 1) and the e- folding number = 60, 
we find ^ = 1674 fixed by the CMB power spectrum 
(procedure explained in the next subsection) . It is found 
numerically that the trajectory is along Imip = 0. We 
thus consider only the real values of if. The infiaton 
trajectory fluctuates in the directions of as shown 
in Fig. [1] where the values of the real and imaginary 
parts of and are plotted along the steepest de- 
scent trajectory. The parameter C for the quartic terms 
in the Kahler potential is chosen to be a moderate value^ 
C = 100 and the fleld values are measured in the reduced 
Planck unit SttG = 1. While the trajectory seems rather 
complicated, inflation actually terminates (that is, one of 



Being higher order terms in the Kahler potential, ^ is expected 
to be not much larger than ^. For zero or small C, the multi- field 
effects become important. While the physics of the isocurvature 
modes and non-Gaussianity arising in the small C, case would also 
be of interest, we shall not discuss such issues in this paper. 



-0.04 




FIG. 1: The steepest descent trajectory of the scalar potential 
([20)). for Ml = M2 = 10" GeV, a = 0, & = 1, 5 = 0, a = and 
= 1, 5 = 1674, = 100. The red solid, the orange dashed, 
the green dot-dashed and the blue dotted curves respectively 
shows the values of ReiVf , ImiVi , ReiV| and ImiVa . The slow 
roll terminates at ^ = 0.0371. The field values are measured 
in the reduced Planck unit SttG = 1. 

Setting Ni, N2 — >■ and restricting (p to take real 
values, the scalar potential (IT71) simplifles to 

Vf = (22) 

where A = {y\jyD)kk is the k-k component {k — 1,2,3) 
of the matrix representation of 

yhvD^ -t^^Umns^'^R^MR^uI^s. (23) 
v' sm p 

The resulting inflationary model is essentially the non- 
minimally coupled X(j)'^ model, discussed in 's^]. We 
have analysed the above example (corresponding to A = 



^ In deriving these values the deviation from N^, N2 = are 
taken into account but the dynamics of A^,"^ are neglected. 
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{y\jyD)ii = 0.00490) in this approximation, and found 
essentially no difference. For example the predictions for 
Us and r are exactly the same within 3 significant digits. 
For smaller Mm the effects of nonzero become even 
more negligible, as the slow roll terminates at a larger 
value of if and the turn around behaviour of occurs 
at the smaller mass scale of Mm- 



2.x 10- 



l.X 10- 




-0.02 



0.04 0.00 



FIG. 2: The scalar potential Ve as a function of Lp and ReTVf . 
The other fields ImA''f , ReA*'!, ImA^'l are chosen to take values 
along the steepest descent trajectory. 



C. Cosmological parameters 

As the single field approximation Nm = 0, € R is 
sufficiently accurate, we shall study the inflationary dy- 
namics within this approximation. We rescale the infla- 
ton field as 



X 



(24) 



so that the kinetic term is canonically normalised in the 
Lagrangian 



1 + ^;^2 I 



(25) 



Here the subscript J indicates the Jordan frame and 



Vf 



16 



4 /■ 

X , ^ 



1 

4 



(26) 



This is the Lagrangian in the Jordan frame where the 
scalar field is nonminimally coupled to the background 
curvature. The Lagrangian in the Einstein frame is ob- 
tained by the Weyl transformation ^Jf^, = (1 -|- ^x^)5^iy 
The canonically normalised inflaton field x in the Ein- 
stein frame is related to x by 



The potential in the Einstein frame is 



A 



X 



16 1 + Cx' 



(28) 



in terms of which the slow roll parameters in the Einstein 
frame are defined as 



1 dVE 
Ve dx 



V = 



1 d^VE 
Ve dx^ ■ 



(29) 



The potential Ve in the single field approximation ([28]) 
is much simpler than the original one (I17p and is deter- 
mined solely by A = {y]jyD)kk and ^. For given A and 
the e- folding number Ng, the nonniinimal coupling ^ is 
fixed by the CMB power spectrum. For definiteness we 
use the maximum likelihood value A\{ko) = 2.43 x 10"^ 
from the 9-year WMAP data 0, 0] , where the calibration 
is set at ko = 0.002 Mpc~^. This is related to the power 
spectrum Vr = VE/247r^e of the curvature perturbation 
(at the horizon exit of the comoving scale) by A|j(fc) = 

^pi'Pnik). We denote the value of x at the end of the 
slow roll (characterised by max(e, rj) = 1) a.s x — X*^ and 
the value of x at the horizon exit of the comoving CMB 
scale A; as X = Xfc- These are related by the e- folding 
number through = J^^' dxVE{dx/dx)/{dVE/dx). 

From the values of the slow roll parameters at the hori- 
zon exit of the comoving CMB scale, the scalar spectral 
index = dlnVR/dhik = 1 — 6e -I- 27; and the tensor- 
to-scalar ratio r = Vgvr/VR = 16e can be computed. The 
results^ are summarised in Table U for — 50, 60 and 
several values of A = {y'^uyojkk- There exists a lower 
bound on A set by the minimal coupling limit ^ — >■ 0. 
In this limit the model reduces to chaotic infiation with 
the quartic potential Ve = ^X^j where A = 1.05 x 10~^^ 
for Ne = 50 and A = 6.19 x IQ-^^ for Ne = 60,^fixed 
by the value of Vr. In contrast to the SM j6l-[l9l| and 
the supersymmetric p7l - l30l | Higgs inflation models the 
nonminimal coupling ^ does not need to be large. We 
see from the table that ^ < for A < 10~^. Hence, 

at least for these values of A the model is obviously free 
from the danger of violating the unitarity bound. Small 
^ is also favoured for avoiding large i?-parity violation, as 
discussed later in Section VI. The Dirac Yukawa coupling 
is related to a, b, 6, a, Mi, M2 through We will 
discuss the relation between A and these parameters in 
Section V. Leptogenesis also constrains certain neutrino 
mass parameters, which will be discussed in Section IV. 



D. Observational constraints and comparison with 
other models based on supersymmetric seesaw 

We show in FigI3]the predicted values of Ug and r of our 
model along with the 68% and 95% confidence level con- 



dx 



yi+^x'+eex' 



dx- 



(27) 



^ These are the results of the tree-level computation. The renor- 
malisation effects are verified to be negligibly small l3Cl |33| . 
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20.2 
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0.165 
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0.106 


0.968 


0.00296 
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0.968 


0.00296 
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0.968 
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0.00303 
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1.18 
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0.968 


0.00369 
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0.00588 
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5 X 10-12 


0.0135 




3.10 


21.5 


0.965 


0.0356 




2 X 10-12 


4.44 X 10" 


-3 


3.32 


22.0 


0.962 


0.0822 




10-12 


1.24 X 10" 


3 


3.42 


22.2 


0.957 


0.161 




6.19 X 10-1^ 







3.46 


22.2 


0.951 


0.260 



TABLE I: The nonminimal coupling ^, the inflaton values at 
the end of the slow roll (x*) and at the horizon exit (xk), the 
spectral index Us, and the tensor-to-scalar ratio r for e- folding 
= 50, 60 and for various different values of A = (ylyyD)kk- 
The coupling ^ is fixed by the amplitude of the curvature 
perturbation. The last hues {Ne = 50, A = 1.05 x 10"i2 and 
Ne = 60, A — 6.19 x 10-1"^) correspond to minimal coupling. 
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FIG. 3: Plot of the scalar spectral index Ua and the tensor-to- 
scalar ratio r, with the 68% and 95% confidence level contours 
from the WMAP9-HeCMB+BAO+i?o data Indicated 
by • are the prediction of our model (NM-I/T^u) for e-folding 
numbers A''e = 50 and 60. The coupling A — {y]jyD)kk is 
a diagonal element of (|23p . which is related to the neutrino 
mass parameters (see Sec. V). For comparison, the spectra 
of the Harrison- Zel'dovich (HZ), the m?<fP minimally coupled 
chaotic inflation model, as well as the A-term-lifed MSSM 
fiat-direction inflation model (AFD) are also indicated. Aws 
is the expected Planck accuracy [5a.[5^. 



[Nf, < 50), we show the case for N^. ~ 50 (thus 
Us = 0.92) in the figure, marked with ^ (AFD) 



tours from the WMAP9-heCMB-HBA0 +Ho data 
The prediction is seen to be consistent with the present 
data, apart from the regions where the Yukawa coupling 
is extremely small. The 2-a constraints roughly corre- 
spond to Mm ^ TeV, details of which depending on the 
other neutrino mass parameters (see Section V). While 
the minimally coupled At/''* model lies outside the 2-(t 
constraints, the nonminimally coupled A^"*, even if the 
coupling is as small as ^ ~ 10~^, sits well inside the 1-a 
constraints. In the figure we also indicate for comparison 
the Harrison-Zel'dovich spectrum as well as the spectra 
of two other inflationary models arising from the super- 
symmetric seesaw model. These are: 

Nr chaotic inflation model — Inflation is driven by a 
right-handed scalar neutrino, with the seesaw mass 
scale identified as the inflaton mass (36l - l38| . It is es- 
sentially the minimally coupled m?(tP' chaotic infla- 
tion model, but has phenomenological advantages 
such as automatic leptogenesis. The spectrum is 
marked with ■ in the figure. 

A-term inflation model — This model assumes 
u'^dPdf-, e'^LL, or N'^LHu direction of the (singlet- 
extended) MSSM as the inflaton [H-Sl]. The 
A-term inflation model predicts very small r and 
~ 1 — 4:/Ne. Since inflation takes place at a 
low energy scale and the e-folding cannot be large 



In the figure we also indicate the expected resolution 
Aus ~ 0.0045 of the Planck satellite experiments (sslls^ l. 
The data from the Planck satellite is soon to be available. 
There are more sensitive CMB polarisation experiments 
planned in the near future (sgI . |57| , which we expect will 
put these models to the test with even higher resolution. 



IV. REHEATING AND BARYOGENESIS 

The scalar potential ([20)1 has a minimum a.t (f = = 
corresponding to the global supersymmetric vacuum. 
After the slow roll, the inflaton undergoes coherent os- 
cillations about this minimum and decays into the SM 
particles, followed by thermalisation. The nonminimal 
coupling can alter the reheating temperature only when 
$ is extremely large and the decay rate is very small 
[ssl . [HI . In the model we are discussing the inflaton is 
in the L-H^ direction and the coupling to the SM par- 
ticles is not small. The nonminimal coupling thus has 
negligible effects on the reheating process. The upper 
bound of the reheating temperature can be estimated as 
Tkh ^ 10'' GeV, assuming the Higgs component decay 



^ There exists a variant of A-term inflation, called inflection point 
inflation [55j While this model allows ris > 0.92, it has less pre- 
dictive power on the spectrum. 
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channel — 56 in the conservative perturbative decay 
scenario. The nonperturbative reheating scenario [60| 
(with parametric resonance) and/or the (s)lepton com- 
ponent decay may lead to a slightly higher upper bound. 
Taking into account the effects of the redshift during the 
time needed for thcrmalisation, we estimate the reheat- 
ing temperature of this scenario to be Trh ~ 10^ — 10^ 
GeV. Below in this section we discuss implication of this 
reheating temperature on the generation of the baryon 
asymmetry. The constraints from the R-parity violation 
will be discussed in Section VI. 

The baryon asymmetry ([1]) can be accounted for by dif- 
ferent mechanisms, depending on the seesaw mass scale 
Mm ■ If the mass scale is smaller than the reheating tem- 
perature Mm ^ Trh, the right-handed (s)neutrinos ther- 
malise and their decay generates lepton number asym- 
metry. The lepton number is later converted into the 
baryon number through the sphaleron process. This 
mechanism is known as thermal leptogenesis 43] . If the 
seesaw mass scale is larger than the reheating tempera- 
ture Mm ^ Trh, on the other hand, the lepton number 
can be generated by the decay of oscillating sneutrinos 
[sgI IstI . l6l| . We call this mechanism nonthermal lepto- 
genesis. The generated lepton number is converted into 
the baryon number, similarly to the thermal leptogene- 
sis case. In addition, as our model includes the MSSM 
components, the AfHeck-Dine mechanism [62] can be op- 
erative. We shall discuss the thermal and the nonthermal 
leptogenesis in our model below. 



A. Thermal leptogenesis 

In the leptogenesis scenario the out-of-equilibrium de- 
cay of the right-handed (s)neutrinos generates the lep- 
ton asymmetry, which is later converted into the baryon 
asymmetry by the [B + L)-violating sphaleron transi- 
tions. In the supersymmetric theory the conversion rate 
is computed to be 



22Nf + ISNh 



23 



(30) 



where Yl is the yield (the ratio of the number density to 
the entropy density) of the leptons, and Nf = 3, Nh = 2 
are the number of the fermion families and the number 
of the Higgs doublets. 

In generic scenarios of leptogenesis, generation of suf- 
ficient baryon asymmetry requires the reheating temper- 
ature to be higher than 10^ GeV. This is much higher 
than the reheating temperature of our scenario, and in 
supersymmetric models the gravitino problem can also 
be serious. It has been pointed out, however, that if the 
two right-handed neutrino masses are nearly degenerate 
the CP-asymmetry parameter is enhanced by resonance 
effects, making the leptogenesis viable even at lower re- 
heating temperature [44 46] . In the following we assume 
the resonant leptogenesis scenario with nearly degenerate 
right-handed neutrino masses Mi w M2. Thermalisation 



of the right-handed neutrinos after reheating also requires 
Ml, M2<Tw 



+ 




+ 




+ 




Hu 



+ 



FIG. 4: The diagrams contributing to the decay mode Nm — > 
The decay modes of the right-handed (s)neutrinos are 



Nr, 



iVil, ^ £,+h, 



(31) 



and the CP-asymmetry parameters associated with these 
processes are defined as 



T{Nm ^ ij +h)- T{Nm -^Il + h) 



r(iV„ ^ +h) + T{Nm -^Il + h) 





T{Nm - 




- T{Nm ' 


^ ^0 + Hi) 




r(iv„, - 


>tj+Hu) 


+ T{Nm - 


-> + Hi) 



11 + Hi) 



r(iV^ ^£,+h)~T{N,-; 



h) 



V- r(7V;^ -^i^+h)+TiN^^£,+h) 



(32) 



Here, denotes the right-handed scalar neutrinos, h 
denotes the (up-type) higssinos, and £, £ are the compo- 
nents of the lepton and scalar lepton doublets. We follow 
the notation of [i^. The CP asymmetry parameters are 
computed from the interference of the tree and one-loop 
diagrams. In Fig. |4]we show the diagrams contributing 

to one of the decay modes Nm —t' £j + h. The resulting 
formula for the CP-asymmetry parameters is 

n^m iyDyh)'mmiyDyD)nn V2 " "/ ' 

(33) 
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where ^M^^ 



Ml 



Ml, and 



M2 



Ml 
MS 



Sfi — 



Ml. AMI 



[AM, 



2 \2 

nmJ ' "'"m-" n 



M2 r2 ' 



(34) 
(35) 



corresponding respectively to the vertex corrections (such 
as the 2nd and 3rd diagrams of Fig. S]) and the self-energy 
corrections (the 4th and 5th diagrams). The decay width 
is written as 



An 



-M„,. 



(36) 



The lepton asymmetry Yl , and hence the baryon asym- 
metry through the conversion pO[) . are found by solving 
the Boltzmann equations. Wc summarise the Boltzmann 
equations and related formulae in the Appendix. Here 
we discuss the results and features that are relevant to 
our model. We first notice that the CP phases 6 and a 
do not affect the baryon asymmetry. This is due to the 
fact that the MNS matrices cancel in the product of the 
Yukawa couplings appearing in (p3| . 



yoVD 



v'^ sin fi 



MRD^RWM. 



(37) 



Here, D^, = diag(m2,m3) for the normal mass hierarchy 
and = diag(mi,m2) for the inverted mass hierarchy. 
The baryon number, on the other hand, depends on the 
parameters a, 6 in (fT3|) . as well as on the right-handed 
neutrino masses Mi, M2. In Figl5]we show the yield of 
the baryon asymmetry Yb , as a and b are varied. The left 
panel shows the result for the normal mass hierarchy with 
Ml = 10^ GeV, and the right panel is for the inverted 
mass hierarchy with Mi = 10^ GeV. In both panels the 
mass difference is chosen to be M2 — Mi = 10~^ x Mi. 
The red contour curves indicate Yb = 8.7 x 10^^^ corre- 
sponding to the observed value. 

The dependance of the baryon asymmetry of the uni- 
verse on these parameters can be understood as follows. 
The baryon asymmetry generated through thermal lep- 
togenesis is proportional to the CP violation parameter 
Em as [64l. 16511 



Yb w k — , 
9* 



(38) 



where sa 200 is the number of degrees of freedom dur- 
ing leptogenesis and k < 1 is the efficiency factor depend- 
ing on details of the Boltzmann equations. In the res- 
onant leptogenesis with nearly degenerate right-handed 
neutrino masses Mi « M2 we find Vn ~ In 2 ^ Sn- 
The maximal enhancement of the CP asymmetry pa- 
rameter Sm occurs when the decay width of either of 
the right-handed (s)neutrinos becomes close to the mass 



difference'^, AM^^ 
difference as AMa 
rameters read 



MrrJOn- Parametrizing the mass 
aMir2, the CP asymmetry pa- 



lm 



iyDyl,)i2 



£2 



1 + "^ {yDyD)ii{yDyl,)22 
(a2 + i)rir2 



-El. 



(39) 



For the (nearly) degenerate right-handed neutrino 
masses, the matrix ^/M in (j37|) becomes proportional 

to the identity. Then Im (yD2/D)?2 /{yDyD)ii{yDyD 
is written as 



22 



2(m 



mg) sin2asinh2& 



(m2 — ^3)2 cos^ 2a — (m2 + TOa)^ cosh^ 2b 
for the normal mass hierarchy and 

2(rnf — m|) sin 2a sinh 2b 

(nil — ?7!,2)2 cos^ 2a — (mi + 7712)^ cosh^ 2b 



(40) 



(41) 



for the inverted mass hierarchy. These have maxima at 
a = 7r/4 w 0.785 and b = ln(3 + 2^2) /i w 0.441, with 
maximal value (ma — 1112) /{m^ + 1112) ~ 0.704 for PO)) 
and (m2 - mi)/(m2 -I- mi) w 7.93 x 10"^ for (gT)). Note 
that Em do not depend on the seesaw mass scale Mm 
except through a, Fi and r2. By adjusting the parame- 
ters a, b and a that are not constrained by observation, 
it is always possible to reproduce the baryon asymmetry 
Yb = 8.7 X 10-". 



-o 1.0 




FIG. 5: The neutrino mass parameters a and b yielding the 
observed value of the baryon asymmetry Yb = 8.7x 10~^^ (the 
red contour curves). The left panel shows the case Mi — 10^ 
GeV with the normal mass hierarchy, and the right panel 
is for Ml = 10^ GeV with the inverted mass hierarchy. In 
both panels the contours are Vs = 2 X 10"^°, 1.5 x 10"^", 
8.7 X 10"", 5 X 10"", 1 X 10"" from inside. We have chosen 
(M2 - Mi)/Mi = 10"^ 



While this might seem enormous fine tuning, it can happen 
naturally as a result of renormalisation group effects [66| |. See 
also il. 
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B. Nonthermal leptogenesis due to decaying 
right-handed scalar neutrinos 



The mechanism described above operates when the 
right-handed neutrino masses Mi, M2 are smaller than 
the reheating temperature which we assume in our sce- 
nario to be in the range Trh = 10^ — 10^ GeV. For larger 
Ml, M2 the thermal leptogenesis scenario is not applica- 
ble as the right-handed (s)neutrinos do not thermalise. 
It is nevertheless possible to generate sufficient baryon 
asymmetry due to the decay of the right-handed sneutri- 
nos that have acquired expectation values during infla- 
tion. In this sense our model is somewhat similar to the 
inflation model driven by the right-handed scalar neu- 
trinos [36l - [38j . Indeed, for large seesaw scales the right- 
handed scalar neutrinos have noticeable contribution to 
the dynamics after the slow-roll, as shown in FiglT] for 
Ml = M2 = 10^^ GeV. Note, however, that the right- 
handed scalar neutrinos do not have to be involved in 
the inflationary dynamics. 

In the supersymmetric SM including the right-handed 
neutrinos leptogenesis is automatic as long as the Hubble 
scale during inflation iJinf is larger than the right-handed 
neutrino mass scale [6l|, which is always the case in our 
model. Take, for example, the case of Mi = M2 = 10^'^ 
GeV and A^e = 60 discussed in Section III B, which gives 
r = 0.00296. Using the CMB power spectrum A|j(fco) = 
2.43 X 10~^, the Hubble parameter during inflation can be 
expressed in terms of the tensor-to-scalar ratio r — IQe 
as ffinf = 2.67 X X lO^'' GeV. With r = 0.00296 
the Hubble parameter during inflation is larger than the 
right-handed neutrino masses, ffinf ~ 1.45 x 10^"^ GeV 
> Ml, M2. The seesaw relation ([SJ with the Yukawa 
coupling < 0(1) constrains Mi, M2 < 10" GeV, and 
for smaller Mi, M2 the model yields larger r and hence 
larger 7?inf, so the inequality iJinf > Mi, M2 is satisfled 
more safely. 

In this nonthermal leptogenesis scenario the right- 
handed scalar neutrinos acquire expectation values {N) , 
either classically by the nontrivial inflaton trajectory, or 
due to the quantum fluctuation during inflation. This 
(iV) can be regarded as the initial value of the right- 
handed scalar neutrinos that oscillate and decay. We 
shall assume {N) < Mpi, where Mpi is the Planck mass. 
Along with the decay of the right-handed scalar neutrinos 
the lepton asymmetry is generated, which is later con- 
verted into the baryon asymmetry by the B — L conserv- 
ing sphaleron process. The process of the decay actually 
depends on the decay rate of the inflaton T^p and that 
of the right-handed neutrinos r„. When (i) the scalar 
neutrino decay rate is larger than the inflaton decay rate 
Tri > r^, the scalar neutrinos decay during the reheat- 
ing. When (h) r„ < < r„(Mpi/(iV))4, the scalar 
neutrinos decay after the reheating but do not dominate 
the energy density of the universe. When (iii) the scalar 
neutrino decay rate is much smaller T^p > r„(A/pi/(A^))'*, 
the right-handed scalar neutrinos dominate the universe 



before they decay. In our model, the scalar neutrino de- 
cay rate^ is (1551) and the inflaton decay rate is estimated 
as « 4 MeV, assuming the Higgs — >■ bb decay channel 
and the Higgs mass ~ 125 GeV. It can be shown that 
in our scenario the case (iii) never occurs, as long as the 
inflaton trajectory is controlled to be nearly straight by 
the noncanonical Kahler terms. The threshold between 
(i) and (ii) is around Mi w M2 ~ 10^'^ GeV, for typical 
values of the mass parameters a = b = 1 and for the 
normal or inverted mass hierarchy. 

In both (i) and (ii) cases the yield of the lepton asym- 
metry is estimated as [6l| 



Yr 



(42) 



and the baryon asymmetry follows from pop. Assuming 



STT^Mi 



(43) 



generated by the quantum fluctuations^, the observed 
baryon asymmetry ([1]) with the condition for the CP 
asymmetry parameter \em\ ^ 1 leads to a mild upper 
bound on the seesaw mass scale Mm 5, 10^^ GeV. 



V. THE INFLATON SELF COUPLING AND 
THE NEUTRINO MASS PARAMETERS 

Our model of inflation is well approximated by the 
nonminimally coupled single fleld A(/)^ model. We have 
seen in Section III that the inflaton self coupling A — 
{y]jyD)kk is related to the CMB spectrum, and argued 
that in the near future the value of A will be constrained 
more severely by the observation. As A is constructed 
from components of the Dirac Yukawa coupling, it is de- 
termined by the neutrino mass parameters. In this sec- 
tion we describe the relation between A and these param- 
eters. 

We shall assume nearly degenerate right-handed neu- 
trino mass parameters Mi « M2, which is favoured for 
successful thermal leptogenesis as we discussed in the pre- 
vious section. The matrix M — diag(Mi,M2) is then 
proportional to the identity matrix and the relation (j23p 
becomes 



yhvD 



2Mi 
v'^ sin^ P 



(44) 



^ Parametric resonance may enhance the decay rate. We do not 

consider such effects here. 
^ If (N) is set by the nontrivial inflaton trajectory as in the case of 

Fig[T] the value of (N) depends on f . Then the lepton number 

generated through coherent oscillation and decay of the scalar 

neutrinos also depends on the parameter 
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Thus {y\)yD)kk is proportional to the seesaw mass pa- 
rameter Ml . Essentially, the inflaton self coupling is de- 
termined by the seesaw mass scale [s^] ■ Further complex- 
ity arises from the dependence on the other parameters 
6, 5 and a, which are not constrained by the present ob- 
servation. 

Dependence on a and b. — We first notice that A does 
not depend on the parameter a. This is easily seen, as 
the product R'^ R appearing in (|44p is written as 

„+ „ / cosh2& isinh26\ 

^'^ =(-zsinh26 cosh 26 j" (^5) 

The couling A monotonically increases with & > 0. We 
show the behaviour of {yljyD)kk/Mi as functions of b in 
Fig. [51 for fc = 1,2,3, with the CP violating phases set 
to he S — 0, a — 0. The self coupling A changes by an 
order as b is shifted by 1; the dependence on b is thus 
significant. 

NH IH 




FIG. 6: Plot of X/Mi = {ylyyD)kk/Mi (in GeV"^) as func- 
tions of b. The Dirac and Majorana phases are chosen to be 
zero. The left panel is for the normal mass hierarchy and the 
right is for the inverted mass hierarchy. On each panel, the 
red solid, the green dashed, and the blue dotted curves indi- 
cate k = 1,2,3, respectively. The k = 2 and fc = 3 curves 
overlie each other. 



Dependence on the Dirac phase S. — The dependence 
of the inflaton self coupling A on the Dirac phase S is 
shown in Fig. [71 where {y]jyD)kk/Mi is plotted for — tt < 
S < TT. We have chosen b — 1 and cr = 0. In the case of 
the normal mass hierarchy the k — 1 component is more 
susceptible than k — 2,3, whereas for the inverted mass 
hierarchy fc = 2, 3 are more susceptible than fc = 1. 

Dependence on the Majorana phase a. — The inflaton 
self coupling A also depends on the Majorana phase a. 
In Fig. [8] we show the behaviour of {y''DyD)kk/Mi as a is 
varied from — tt to tt. We have chosen 6=1 and 6 = 0. 

To summarise, the inflaton self coupling A depends on 
b, S and a but not on a. These parameters are not 
constrained by present experiments and thus introduce 
ambiguity in prediction of the model. Consider for in- 
stance the 2-(T contour in Fig[3l For = 60 this gives a 
constraint A = {y£)yD)kk ^ 10^^^ which corresponds to 
Ml « Af2 > 2 TeV for 6 = 1, 5 = cr = 0, fc = 1 in the 
normal mass hierarchy (Fig[7]). If we take S = 1.6 instead 
of (5 = 0, the 2-cr of CMB gives Mi w M2 > 20 TeV. 



NH 




FIG. 7: Plot of A/Ml = {ylyyD)kk/Mi for fc = 1 (the red 
solid curves), k = 2 (green dashed) and A: = 3 (blue dotted), as 
functions of the Dirac phase 5. The unit of A/Mi is in GeV~^. 
The left panel shows the case for the normal mass hierarchy 
and the right panel is for the inverted mass hierarchy. We 
have chosen 6 = 1 and u = 0. 




FIG. 8; Plot of A/Ml = iyl}yD)kk/Mi for fc = 1 (the red 
solid curves), k = 2 (green dashed), fc = 3 (blue dotted), 
as functions of the Majorana phase a. The unit of \/M\ 
is in GeV~^. The left panel shows the case for the normal 
mass hierarchy and the right panel is for the inverted mass 
hierarchy. We have chosen 6=1 and (5 = 0. 



VI. THE EFFECTS OF R-PARITY VIOLATION 

While the seesaw model superpotential ^ preserves 
the R-parity (assuming the odd parity for the right- 
handed neutrino superfield), the Kahler potential (fT6)) 
breaks it by the "/iLiHu terms. Small R-parity viola- 
tion is known to be harmless, and it can often be bene- 
ficial [68l - l70j . It can nevertheless lead to difficulties and 
the consequences of introducing such terms need to be 
checked carefully. In this section we discuss viability of 
the scenario focusing on the effects of R-parity violation. 
We assume the canonical form of the Kahler terms for all 
the other chiral multiplets. 

An immediate consequence of the 7iLiiJ„ terms in the 
Kahler potential is that the the Lagrangian includes the 
following terms: 

C D J dW4'{LlL, + HlHu+j.,L,H.a + h.c}j 

= J d^d + HlHu + Jl^^H,, + h.c}j (46) 

where (j) is the compensator in the superconformal for- 
malism. To go from the first to the second line we 
have rescaled Li,H^ — >• (j)~^Li,(l)~^H^. When the su- 
persymmetry is broken the compensator acquires an ex- 
pectation value (0) =1-1- O'^Fij,, where is the com- 
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pensator F-term. The third term in (|46)) then becomes 

/ d^eil + 5VJ)(1 - e^F^)j,UH^. This suggests gen- 
eration of a bihnear R-parity breaking term Wfi = 
Pj^'jiLiHu = fiiLiHa, as well as an R-parity breaking 

B-term —\F^\'^ji£ihu. In generic scenarios of supersym- 
metry breaking® the compensator F-term gives rise to the 
gravitino mass, « Recalling that 7^ = 7 is re- 

lated to the nonminimal coupling C by ^ = ^ — ^, we see 
7 w 45 > 1 when C > 1 and 7 0(1) when ^ < 1. We 
will assume 7 ^ 0{1) below, as one of the merits of our 
model was that the extremely large nonminimal coupling 
that could lead to the unitarity violation can be avoided. 
The bilinear R-parity breaking term can then be written 
as Wft ~ m3/2jiL.iHu m^i2LkHu. 

In the presence of the bilinear R-parity violating terms 
W ^ Hu{^J.Hd + m^/2Lk) the lepton number violating 
superpotential 

I^AL=i ~ {yl'e'')e'rL,Lk + {y^^ e'')d1Q,Lk, (47) 

where ^ "^3/2/ is generated. The cosmological con- 
straints for these effective Yukawa couplings are (72| - [7^ 

y'^e^ 2/y£'=<10-^ (48) 

from which one obtains ^ ™3/2/m ^ 10^^ cos /3 ^ 
10^^. For a typical value of the MSSM ^ parameter 
/i ~ 1 TeV the gravitino mass of our model is 

^3/2 < 1 MeV. (49) 

Thus, supersymmetry breaking mechanisms consistent 
with our scenario are those giving a small gravitino mass 
(such as the gauge mediation scenario). Note that the 
gravitinos with the small mass (|49p are still a good can- 
didate for the cold dark matter, even though the R-parity 
is broken. 

The abundance of the thermally produced gravitinos 
is calculated as [75l-l77| 

^ V 10^" GeV J V m3/2 J \l TeV J ' 

(50) 

where h is the Hubble parameter in the unit of 100 km 
Mpc~^ s~^ and M-g is the running gluino mass. With 
W3/2 ~ 1 MeV and Mg « 1 TeV, the condition for 
avoiding the overdominance fl^h^ ^0.1 yields the upper 
bound of the reheating temperature Trh ^ 10^ GeV. In a 
previous section we estimated the reheating temperature 
be Trh ~ 10^ — 10^ GeV, assuming the decay channel 
Higgs — bb. The constraint of the gravitino abundance 
(|50)) suggests that Trh ~ 10^ GeV is favoured. Taking 
into account the Hubble expansion for the time needed 



There are exceptional cases, e.g. the almost no-scale model |71| . 



for thermalisation, Trh '-^ 10^ GeV seems to be a rea- 
sonable reheating temperature of our scenario. To sum- 
marise, our model typically predicts a scenario of grav- 
itino cold dark matter with the reheating temperature 
Trh ^ 10^ GeV. The smallness of the R-parity violat- 
ing terms (j48|) guarantees that the baryon asymmetry 
generated by thermal/nonthermal leptogenesis remains 
without being wiped out. 

We close this section with two comments. First, the 
constraint on the R-parity violation does not have 
to be taken too strict since, for example, details of the 
flavour structure may slightly relax this condition [ts} . 
Our second comment concerns the effects on the neu- 
trino masses. R-parity violation induces neutrino masses 
without invoking right-handed neutrinos [79l - [83 |. This 
alternative mechanism to the seesaw can be a poten- 
tial threat, since if such an effect dominates over the 
usual seesaw mechanism that would give unacceptably 
large neutrino masses. A back-of-the-envelope calcula- 
tion shows that there is actually no such danger, as the 
constraint (|48l) is more stringent than the one coming 
from the neutrino masses. For example, in the bilinear 
R-parity breaking scenario [sH - fSTf the condition on the 
bilinear coefficient /x^ for not generating too large neu- 
trino masses is Mi/yU < 10"'^ >• 10^^. 



VII. DISCUSSION 

We have discussed in this paper a simple model of infla- 
tionary cosmology based on the supergravity-embedded 
minimal seesaw model. The scenario is economical as it 
simultaneously explains various issues - neutrino oscilla- 
tion, the origin of the baryon asymmetry and the origin 
of the cold dark matter ~ apart from the standard is- 
sues of big bang cosmology solved by inflation. We have 
shown that the model reproduces observationally accept- 
able values of the cosmological parameters, and argued 
that the prediction for rig, r of the CMB spectrum will 
be tested by satellite experiments in the near future. A 
particularly interesting feature of this model is that the 
seesaw mass scale is constrained by the CMB. Thus far 
useful constraints on the (left-handed) neutrinos, such as 
the total neutrino masses ^ nii, and the effective num- 
ber of the neutrino species, have been obtained by ob- 
serving the CMB and the large scale structure. However, 
the nature of the right-handed neutrons remains myste- 
rious: being gauge singlets, their detection in colliders 
is virtually impossible, nevertheless they are essential for 
both seesaw mechanism and leptogenesis. In our pro- 
posal, the CMB may provide access to the physics of the 
right-handed neutrinos. 

A key feature of our model is the nonminimal coupling 
of the D-flat direction inflaton to the background grav- 
itational curvature, which is naturally implemented by 
supergravity embedding of the SM. In contrast to the 
nonmimimally coupled Higgs inflation type models, the 
coupling in our case need not be large. This is related 
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to the fact that the Dirac Yukawa couphng can be very 
small. Such extremely small Yukawa coupling is, never- 
theless, not unnatural. The Dirac Yukawa coupling cor- 
responding to a TeV scale right-handed neutrino mass in 
our model is in the same order as the electron Yukawa 
coupling. As Nature allows such a small coupling for the 
electrons, there is no reason to exclude the Dirac Yukawa 
coupling of the same order for the neutrinos. 

Finally, we comment on extension of our model in var- 
ious directions. While embedding into the 50(10) grand 
unified theory is not possible, one may for example con- 
sider our scenario in the grand unified theory of SU (5) 
plus a singlet. Also, type III seesaw with SU{5) adjoint 
neutrinos in SU{5) is possible. It is also straightforward 
to extend our model to the right-handed neutrinos with 
three families. An obvious drawback of such an exten- 
sion is that the inflationary scenario will contain more 
unconstrained parameters and the predictive power of 
the model will be reduced. 
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Appendix A: Boltzmann equations 



lepton: 



H,, 



h + e. 



The decay widths for these processes are, 



47r 



iyDyD)r 



(A2) 



(A3) 



It is convenient to parametrize the inverse temperature 
using the right-handed neutrino mass Mi as 



Ml 



(A4) 



In the out-of-equilibrium decay of the right-handed neu- 
trinos and sneutrinos, we assume the initial number den- 
sities of the leptons and sleptons Yl^ to be zero. 
We also assume that the right-handed sneutrinos are ini- 
tially symmetric: = Y~-t. Ignoring subdominant 

processes, it follows that = i^^t persists during the 
subsequent evolution. The Boltzmann equations for the 
number densities of the right-handed (s)neutrinos and 
the (s)lepton numbers then read 



dYN„ 
dz 



dY; 



dz 
dz 



dYL^ 
dz 



z fYN„^ \ 

sHiMi) [Yg j^^-' 



(A5) 
(A6) 



1N„ 



], (A7) 



sH{M_ 



1N„ 



1 / ^N. 



2 I Y^'i ' r^'^ 



]■ (A8) 



We discussed leptogenesis within our inflationary sce- 
nario in Sec. IIV Al and presented the solutions of the 
Boltzmann equations. In this appendix we collect re- 
lated formulae. The Boltzmann equations in the context 
of leptogenesis are discussed in j88l - l90| . We follow the 
conventions of [6^. 

In the supersymmetric minimal seesaw model the dom- 
inant processes for generating the lepton number are the 
decay of the right-handed Majorana neutrino N„i into the 
up-type Higgs boson and a lepton, or into the higgsino 
and a scalar lepton: 



Here, H{Mi) is the Hubble parameter at temperature 

T = Ml and 



1N„ 



eq Kl{z) 

N„ ~- ' " 



K2{Z) 



(A9) 



is the reaction density of the decay processes. Ki, Ki 
are the elliptic integrals of the first and the second kind. 
The Bohzmann equations (IX5|) and ([M)) . (|T7t and 
are identical due to supersymmetry. 

We solved the above equations with the yields in equi- 
librium. 



N^^Hu+e, h + i, 


(Al) 




_ ycq 


cq 
S 


as well as the decay of its superpartner 
gsino and a lepton, or into the Higgs boson 


into the hig- 
and a scalar 




e 


cq 

s 



cq _ 



M? 

^K2{z), (AlO) 
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These Boltzmann equations assume supersymmetry and symmetric case is however expected to be minor as it 
are not strictly apphcable below the supersymmetry should naturally be within a factor of 2. 
breaking scale T ~ TeV. The deviation from the super- 
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